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1-009 ° ” 9 on oooh Oke 01:36:31 EATHRTL . SREIATHTANH OMAR: 1 ° (1) 
-TITLE MTHSTANH : F heating Point Hyperbolic Tangent routine 


IDENT /1-009/ ; File: MTHTANH.MAR Edit: JCW1009 
PITPIS ISP iii iii iii ii iiiiiiiiiiiiiiiiiiiiiititiiiiiitiiiiiitititt) | 
COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 

RIGHTS RESERVED. 


FTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIE 
NCE WIT E F SUCH LICENSE AND WITH TH 
N OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
HEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
RSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
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FACILITY: MATH LIBRARY 

+ 

ABSTRACT: 

MTHS$STANH is a function which returns the floating point prperkel ve sengent 


of its single precision floating point argument. The call is standar 
call-by-reference. 


VERSION: 01 
HISTORY: 
AUTHOR: 
Peter Yuo, 29-Jun-77: Version 01 


MODIFIED BY: 
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1-009 HISTORY T petailed. Current Edit History ~$Fb=1 38 94:35:34 MTHRTL.SR CIM THTANH.MAR; 1 3) | 
9 -SBTTL HISTORY ; Detailed Current Edit History | 
; ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none | 
8 5 ; Edit History for Version 01 of MTHS$TANH 
0 5 3 o-4 - Call MTHSCOSH, MTHSSINH directly spetend ¢¢ of using POLY. TNH 13-June-78 
0 8; 0-5 = Declare xternals, word offsets. TNH 13-June-7 
000 59 ; 1-006 = Put seve on number in standard format, with three, di its o 
000 60 ; edit — and ePair copyrigns. nosice. n° novira” 
900 : : 1-007 = Add 5 aot the P PSECT d re tive. 
; 1-008 - Change snet ont .0 oe r, improved. sii icloaey, The 
0000 635; alue of X above oh th 1° 0 i est machine approximation 
0000 64 ; to TANH(X) is about 9.01. the’ next higher number that can be 
0000 65 ; represented as a short literal is 10.0. JAW 19-Sep-80 
0000 66 ; 1-009 - onan ed the constant 2°-14 to 2%-12 te correct inaccuracy. For 
0000 67 ; ues of {Xi between Deals and 2*-12 no noticeable loss of 
0000 68 ; significance, wes noticed by the assumption that TANH(x)=x starting 
0000 69 ; at iXi<=2*-14 instead of 2 but a loss of performance was felt 
0000 70 ; do to an unnecessary cacoueee ibe of SINH(X)/COSH(X). ALL appropriate 
44 | : references to 2“-14 have been changed to 2“-12. JCW 10-Jan=83 
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1-009 BECLARATIONS. Declarative Part of Modul g-SEp 7138 94:39:35 ERTMRTL. SREIMTHTANH OMAR: 1 ei | 
73 -SBTTL DECLARATIONS : Declarative Part of Module 
r 3; INCLUDE FILES: OTSPARAMS .MAR 
5 3; EXTERNAL SYMBOLS: MTHSJACKET_HDLR 
Be .DSA ABL GBL : Force all symbols to need .EXTERN 
3 EXTRN nTisexP RG : EXP " | 
0 “EXTRN MTHSCOSA ; COSH 
00 4 »EXTRN TSS aN 3; SINH | 
$8 ‘ ; EQUATED SYMBOLS: | 
0 00 .25 = 0816 + _“X3F 80 : 0.25 
80 PAB 00 8 tf -fOWER,, M1 = 0a16 + *X3A80 ; 952240 
990040 4 3 9 pa *F1.0 : 1.0 
00004 0 90 5F710.0 = *F10.0 ;.10. 
00000004 0 91 value =4 : value.rf.r 
4 3 
000 93 ; 
0000 94 ; MACROS: none 
44 95 ; 
000 38 3; PSECT DECLARATIONS: 
0000 g 
00000000 98 -PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
0000 99 3 program section for math routines | 
0000 100 ; 
0000 101 ; OWN STORAGE: none 
0000 19 3 
0000 103; 
9000 104 ; CONSTANTS: 
0000 105; 
0000 106 ; NONE 
| 
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Nn Fl 6=SEP=1984 MTHRTL.SRCIMTHTANH.MAR; 1 
-SBTTL MTHSTANH = Standard Single Precision Floating TANH 


FUNCTIONAL DESCRIPTION: 
TANH = single precision floating point function 
TANH(X) is computed as: 

If iX! =< 2ee-12, then TANH(X) = X, 


If gesct2 < 1X! =< 0.25, then TANH(X) = SINH(X)/COSH(X). 
If 0.25 < 1x} < 10.0. then TANH(X) = (EXP(2*x) = 1) / (EXP(2*x) + 1) 
If 10.0 =< Xi, then TANH(X) = sign(X) * 1 


CALLING SEQUENCE: 
TANH.wf.v = MTHSTANH(x.rf.r) 
INPUT PARAMETERS: 


LONG = 4 ; define eepere multiplier 
x = 1 * LONG ; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: floating hyperbolic tangent of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: none 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


«ENTRY MTHSTANH, “M<IV, R2, R3, R4> 

; standard call-by-reference entry 
disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


MOVAB G“MTHSSJACKET_HND, (FP) 


set handler address to jacket 
handler 


; case of an error in routine 
; If an error, convert signal to user PC 
3 and resigna 


MOVE @value(AP), RO Q = IX! = @value(AP) 
BiCw  #*x8000, RO RO = iX! 
CMPF RO, #LF_2_POWER_M12 compare Xi with 2*-12 


dD 3 
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1-009 Hs TAN 9 Standard Single Prevision ft raat 94:38:34 YRTARTE SREDATHTANN MAR: 1 — (ay I= 
46 «15 1? 199 BLEQ OUT_X ; branch if IX! =< 2ee-12 
18 196 3 
18 «165 5 2e*-12 < IK! 
18 164; 
22 51 ig 182 CMPF RO, S*#SF_10.0 x ith 10.0 
“#SF 10. ; compare {Xi wit . 
3 18 5 187 BGEQ GEG_T0_1070 i branch ff ix: >= 10.0 | 
9 169 : 
170 ; 2e*-12 < ix! < 10.0 
6 173 ; | 
Q00003F80 BF 5 51 00 17 CMPF RO, #LF_0.25 z compare !X! with 0.25 
39 15 00 : 174 BLEQ e610, 0038 i branch if ixt =< 0.95 | 
0029 175 | 
0029 176 3 
0029 177 ; 0.25 < ix! < 10.0 
0029 178; 
0029 179 
50 04 BC 04 BC 41 bo5e 139 ADDF3 a@value(AP), avalue(AP), a oer | 
+4 = ® | 
00000000'GF 16 O03 182 SB G*MTHSEXP_RG ; RO = EXP(2*x) 
51 50 08 41 0035 18 ADDF3 S*#SF_1.07 RO, RI ; RI = EXP (geX) +1 
50 08 42 0039 184 UBF S*#SF71.0, RO = EXP(S#x) = 1 | 
50. 5 46 003C 185 DIVF » RO 3 RO = (EXP(2*xX) = 1) / CEXP(2*x) + 1) 
04 OO3F 186 ET 3; return with result in RO 
0040 187 
0040 188; 
0040 189 ; 2-12 < IK! =< 0.25 
0040 190; 
0040 191 | 
0040 13 LEQ_T0_0.25: 
00000000 ' GF 6C FA 74 19 CALLG (AP), G*MTHSCOSH ; RO = COSH(X) 
2 50 50 004 194 MOVF ‘ 3 R2 = COSH(X) 
00000000 * GF 6C FA OO4A 195 CALLG (AP), G*MTHSSINH ; RO = SINH(X) 
50 52 46 0051 196 DIVF  R2, RO > RO = SINH(X) 7 COSH(X) 
04 0054 #4197 RET : return with result in RO 
0055 198 
0055 199; 
055 00 ; ix! >= 10.0 
055 01; 
055 8 
B22 GEQ_T0_10.0: 
50 0O 29 5 4 MOVF S“*#SF_1.0, RO ; RO = 1.0 
048C 53 0058 205 TSTF  a@value APs : test the sign of x 
18 058 Be BGEO 10$ ; branch if xX >= 0 
50 5 05D 0 MNEGF RO, RO 3 =- 
04 0060 88 10$: RET 3; return with result in RO 
061 9 
sd 10 ; 
61 11 5 Iki =< 2ee-12 
as 
50 04 BC 50 61 14 OUT_X MOVF @value(AP), RO ; RO = TANH(X) = X | 
04 2 1? R : return with result in RO 
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Symbol table se ° ut eee a1b4 $1:59:35 FATANTE eReoRritAmeman:1 | 000 
GEQ_ 10.10 SSR 1 
LEG"T -9°38 40k 8} 
LF. 5 = F 
LF 2 _POWER_Mi2 = A 
LONG = 
MTHSSJACKET_HND teeeenee x 6} 
MTHSCOSH eeeeeeee YX 9 
MTHSEXP_R4 renteree x 
MTHSSINA rereeene x 09 
MTHSTANH 8 0000 RG 1 
QUT x BY 1R 01 
= 

5F710.0 = Gug42 8 
vACUE = 00000004 

tawmeose ner ecomocae + 

! Psect synopsis ! 

geeoeeeeecoeeceoceooecs + 
PSECT name Allocation PSECT No. Attributes 
- «ABS. 00000000 <( 0.) OO ¢ QO.) NOPIC USR CON ABS ~ LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 00000066 ¢ 102.) O1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 

Swe ewe new ene ee weer ere wmoece & 
: Performance indicators H 

Phase Page faults CPU Time Elapsed Time 
Initialization 29 80:80:00 .18 00:00:01.48 
pears erouneee "bg HE88 8898 HE a8 96 
éss : : ° : : ° 
pyabo} table sort 53 0:00: 8-92 0:00;08-08 
ass :00:00. :00:03. 
Symbol table output 3 8: :00.0 B3 #20: 0.0 
Psect synopsis output B38: 0.0 B09: 88-8 
Cross-reference output 1: 3-9 8 :00:00.0 
Assembler run totals 294 0:00:02.1 0:00:12.6 


3329 bytes (7 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 15 non-local and 1 local symbols. 
79 source Lines were read in Pass 1, producing 11 object records in Pass 2. 
page uf virtual memory was used to define 1 macro. 


The working get Limit was 900 pages. 
0 


+ 
! Macro Library statistics H 
deem wenn wer ene rene momen oe 


Macro library name Macros defined 
-$255$DUA28: CSYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


ee " 


MTHS 


TANH MT 
VAK-11 Macro Run Statistics 
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MACRO/ENABLE=SUPPRESSION/D] SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHTANH/OBJ=OBJ$:MTHTANH MSRCS$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) #MSRCS: 
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